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Abstract 

In this talk we present techniques for calculating one-loop amplitudes for 
multi-leg processes using Feynman diagrammatic methods in a semi-algebraic 
context. Our approach combines the advantages of the different methods allow- 
ing for a fast evaluation of the amplitude while monitoring the numerical stability 
of the calculation. In phase space regions close to singular kinematics we use a 
method avoiding spurious Gram determinants in the calculation. As an applica- 
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tion of our approach we report on the status of the calculation of the amplitude 
for the process pp — > bbbb + X. 

1 Introduction 

The experimental precision expected to be reached at the LHC demands cross-section 
calculations to move from the leading order approximation (LO) to the next-to-leading 
order (NLO) for many processes with three, four or even more particles in the final 
state. The complexity introduced by both the real and virtual corrections can only be 
overcome using automatised tools in the construction of programs for the numerical 
evaluation of their matrix elements. 

An important step towards the automatisation of NLO calculations was the formula- 
tion of local subtraction terms which render real and virtual corrections independently 
infrared finite. The two most commonly used subtractions at the one-loop level are 
the dipole formalism [TJ[2] and the FKS subtraction method [3JH], both of which have 
led to several packages for their automated implementation [5HT2] . As one of their ad- 
vantages, subtraction methods offer a well-defined way of separating the different parts 
of a one-loop calculation. Inspired by this separation, at the Les Houches workshop a 
standard interface between Monte Carlo tools, computing the Born level, the real cor- 
rections and the subtraction terms, and one-loop programs which only need to provide 
the virtual corrections has been proposed [T3] . 

The persistent efforts spent on the improvement of the methods for the calculation 
of one-loop amplitudes have recently induced a rapid growth in the number of NLO 
predictions for the LHC. The most recent results with four and more final state particles 
at the LHC are the corrections to pp — > bbtt + X [TH4T6] . pp — > V + 3j [TTHT9] . pp — )■ 
W ± + Aj |20], pp -> W + W + jj [2TJ and the qq initiated channel of pp -> bbbb + X [22]. 
However, many Standard Model processes of similar complexity which are relevant for 
the discovery and/or the measurement of the properties of the Higgs boson or new 
particles from extensions beyond the Standard Model (BSM) are not yet available at 
NLO [23]. Moreover, once new particles are found at the LHC, their detailed study will 
require NLO studies in BSM models. We have therefore developed and implemented 
methods for the automated evaluation of one-loop Feynman diagrams, which we de- 
scribe in Section |2j In Section [31 we present applications of our method, one of which 
is the calculation of the gluon induced channel of the process pp — > bbbb + X. 
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2 Method 



2.1 A Numerator Representation for Feynman Diagrams 

Our method for the calculation of one-loop matrix elements is based on the evaluation 
of Feynman diagrams. A one-loop diagram can be written as 

fd»q A%£) 

J in»/ 2 Dx---D N ' { ) 

where Di = [(q + rj) 2 — mf + i8] and the integration momentum can be split into its 
projection <f M on the physical Minkowski space and an orthogonal component q. Instead 
of a single diagram, T> can also refer to a group of diagrams sharing the same set of 
denominators. All quantities are defined such that the following equations hold: 

9^ = 9^ + 9^, ^ = 4, ££ = -2e, (2a) 
g»9 pv = r\ 9$9 pv = r', 9^ = 0, (2b) 
<f = g»q\ ~q» = ~g»q\ q 2 = Q 2 -^- (2c) 

The numerator of the integral can then be expressed in terms of the quantities (q, /i 2 ) 
as 

N{q, e) = N Q {q, /i 2 ) + eMi{q, li 2 ) + e 2 Af 2 (q, fi 2 ) + 0(e 3 ). (3) 

The functions H\ and M2 are relevant in regularisation schemes, where internal gauge 
bosons are kept in n dimensions (n — 2 polarisations), such as the 't Hooft-Veltman 
scheme or naive dimensional regularisation. The terms originating from Af 2 are of purely 
infrared origin and therefore have to cancel out, as shown in [13]. These terms can serve 
as an additional check on the amplitude or simply be skipped in the computation. In 
the following discussion there is no need to distinguish between the individual functions 
jVi(q, /i 2 ), as we can simply consider the calculation of a diagram as separate calculations 
of the three terms in T> = V Q + eD\ + e 2 V 2 . We refer to these functions collectively 
as Af(q, fi 2 ). 

For the algebraic construction of M{q, n 2 ) from the underlying Feynman rules we use 
QGraf (23] for generating the diagrams and Form [25] for their algebraic simplification. 
Since we work with helicity projections of amplitudes, we have devised a Form library 
for the manipulation of expressions containing helicity spinors [26] which is publicly 
available^. 

In amplitudes containing external partons, the matrix element and therefore each 
diagram is a vector in colour space T> = ^2 C T> C |c). We avoid working with the vector 
components T> c by contracting each diagram with the tree level amplitude, and working 
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with the corresponding numerator functions 

A/f CD) (g,/i 2 ) = £ (AL^ K C (Q^ 2 ) (c\c) • (4) 

c,c' 

This step is important to avoid a proliferation of numerical tensor reductions. In fact, 
one cannot avoid working with the components T> c for processes where the Born level 
vanishes and at the first non- vanishing order one has to calculate the square of the NLO 
matrix element!!. 

2.2 Tensor Reduction at the Integrand Level 

Complete reducibility of one-loop diagrams to scalar integrals implies that the g-dependence 
of the numerator can be expressed in terms of the denominators Di 

Ntf, l?) = a h-3N-i (q) D 3i ■ ■ ■ D 3N-x 
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The calculation of a one-loop diagram therefore amounts to the extraction of the co- 
efficients, which are the residues of the integrand at the cuts Dj 1 = . . . = Dj r = 0. A 
suitable parametrisation of the coefficients and a fast implementation of the method in 
terms of a Fortran 90 library S AMURaJI is described in [32] • 

The decomposition given in Equation ([5]) works with n-dimensional rather than 
four-dimensional propagators Di = [(q + rj) 2 — mf\ and keeps the dependence on /x 2 
in the numerator. This improvement over the original OPP method [27J[28] allows the 
determination of the amplitude including both the cut-constructible and the rational 
part. The Samurai implementation also includes several reconstruction tests which 
provide a measure for the quality of the numerical result for a given diagram. This in- 
formation can be used to dynamically switch to an alternative reduction method when 



2 We keep the name NLO for the one- loop order even if the Born level vanishes and the one- loop 
matrix element formally becomes the LO. 
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the precision of the reconstruction is not sufficient. Typically, unitarity based reduction 
methods fail in the neighbourhood of kinematical points with vanishing Gram determi- 
nants. Traditionally, implementations of unitarity methods go to higher precision when 
such a case is detected. The following two sections introduce a method which avoids 
the use of higher precision implementations. 

2.3 Improved Algebraic/Numerical Tensor Reduction 

The shortcomings of tensor reduction methods in the neighbourhood of vanishing Gram 
determinants can be cured either by expanding tensor coefficients about limits of van- 
ishing Gram determinants [33J or by extending the integral basis by integrals which are 
non-scalar but which have fast and stable implementations [33H36] . The latter approach 
has been followed in implementing the Fortran 95 library Golem95 |37fi which in its 
current version can be used for the reduction of tensor integrals with up to six exter- 
nal legs for massive and massless loop-propagators. For some R < N, the numerator 
structure of a Feynman diagram in a gauge theory can always be written as 

/Jn S^R /^( r ) n fj,i . . . n Hr JiL 
zW2 D l ---D N -Z^i.-./v'iV , W 

r=0 

where the tensor integrals have the decomposition 

T n,m...iM- _ ST W- ■ ■ -r- 1^-^ A N ' r + V" \n"r- ■■■r- 1 W '" Mr R N ' r 

J iV - ^ ^\ A h-3r + 2^ I 9 T 3l r >"2J ''< -J, , 

OX— Or OX—Or-l 

+ E \9--9-r n ---r Jr _r^C^, (7) 

JX—jr-4 

The coefficient functions A ,r , B ,r and C N,r are implemented in the Golem95 library, 
which apart from these tensor coefficients, also provides all the scalar integrals, both 
finite and divergent one^], involving massless and massive particles with real masses. 
An extension of the library which, amongst other improvements, allows for complex 
masses will be released soon. 

2.4 Tensorial Reconstruction at the Integrand Level 

In Equations (jSJ) and ([6]) we introduced two different decompositions of the numerator 
function Af(q,fi 2 ), which can be rewritten as 

L-R/2J R-2a 

a=Q r=0 

4 http : / /lappth. in2p3 . f r/Golem/golem95 .html 

5 regulated dimensionally 
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It turns out that a method for the numerical reconstruction of the tensor coefficients 

C (r,a) 

can be constructed by means of extracting the coefficients of a multivariate poly- 
nomial in the variables (q°, . . . , g 3 , /z 2 ). As explained in Section [2^2| a numerical imple- 
mentation of J\f(q, fj 2 ) is sufficient as input for such a reconstruction, which is described 
in detail in [38J . For nearly all relevant cases tensorial reconstruction (jBJ) requires fewer 
samplings of the numerator function than reduction at the integrand level (jSJ). An im- 
provement of the Samurai algorithm can therefore be achieved by passing a numerator 
function to the tensor reduction which has been obtained from the original numerator 
by tensorial reconstruction. Another possibility is to combine Samurai and Golem95 
as shown in Figure [U 



Construct Tensor Coefficients 






I 




Run Samurai with reconstructed N(q) 






ailure 






^^^i^econsiiucuon rest 
success! 


Contract Tensor 


Coeff. in golem95 



Figure 1: Possible ways of using Samurai and Golem95 together with tensorial re- 
construction. On the left branch Samurai is called with the numerator reconstructed 
according to Eq. OH]), which can speed up the calculation. If the consistency check in 
Samurai fails Golem95 is called as a rescue system. On the right branch the call to 
Samurai is skipped and Golem95 is used as the sole tensor reduction method. 



3 Applications 

3.1 Automated One-Loop Matrix Element Evaluation 

We have combined the components described in Section [2] to build an automated gen- 
eral one-loop evaluator for matrix elements (GOLEM). The implementation is realised 
in form of a python program called golem-2.0 [39T44T] . We have recently added an 
interface for importing Feynman rules from FeynRules [12]. Furthermore, the program 
supports the Binoth Les Houches Accord [13] which defines a standard way of com- 
municating information between Monte Carlo generators and one-loop matrix element 
codes. The numerical code is generated in Fortran 95 using the optimising code gen- 
erator haggies [13] but can be linked both by the Fortran [H] and C conventions. 
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3.2 Results for pp ->• bbbb + X 



At the LHC the signature of four 6-jets can be a significant background for Higgs 
searches in models beyond the Standard Model [231H51H6] . Owing to large uncertainties 
and huge backgrounds this search channel has been neglected in most phenomenological 
studies for the LHC. Since in some BSM scenarios one of the Higgs bosons can be 
predominantly decay into a four-6 final state, it will be crucial to include the NLO 
corrections in a s in future studies. 

In our setup, the virtual matrix elements have been generated using golem-2 . for 
the processes qq/gg — > <?i<?i<?2<72, where q, q± and q<i represent different quark flavours. 
The amplitudes for the bbbb final state have been obtained using the relation (40] 

A bM {l, 2; 3, 4, 5, 6) = A qi g iq2ij2 (l, 2; 3, 4, 5, 6) - A^^l, 2; 3, 6, 5, 4), (9) 

where a permutation of the momenta encodes a simultaneous permutation of momenta, 
colour and helicity labels. This trick reduces the number of diagrams and therefore the 
code size by a factor two. The Born level amplitude and the real emission corrections 
were calculated using MadGraph [47] and MadEvent [IS]; for the subtraction of the 
infrared singularities MadDipole [HUTU] has been employed. On top of the programs 
described in Section |2j the library OneLOop [49] was used for the evaluation of the scalar 
integrals. The virtual corrections were integrated over phase space by reweighting sets 
of unweighted Born level events {pi}- In order to obtain the virtual corrections for an 
observable defined by the measurement function ^{{pt}) we evaluate the Monte Carlo 
sum 

(^virt = ^ £ ( 1 + A *™ A *« + k C ) T({Pi})- (10) 
t=l \ ^Born^Born / 

In order to obtain an estimate for the computational cost and the numerical stability 
of the generated code, we have evaluated the matrix element for 500 points, starting 
from a fixed kinematics as specified in J5U] and rotating the final state particles about 
the z-axis which coincides with the direction of the incoming particles. 

The results for the gluon induced subprocess are shown in Figure [2j As expected [5QI 
I5T] the amplitude develops a relatively sharp peak around 9o ~ 2.32 and at Q\ = 9 + n 
because these points are close to a double parton scattering kinematics. On a Xeon 5500 
CPU the computation of a single phase space point in our setup requires 18 seconds for 
the virtual corrections and 11 ms for the real corrections including subtractions. 

4 Conclusion 

In this talk we have reported on the status and various developments in the Golem 
and the Samurai project. Samurai is a program for the unitarity based reduction of 
one-loop amplitudes which, amongst other improvements, preserves the n-dimensional 
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rotation angle 9 



Figure 2: Born (dashed) matrix element squared and the interference term between 
Born and virtual corrections (straight line) of the gg — > bbbb amplitude. The curves 
have been obtained by rotating the final state particles about the z-axis. The amplitude 
peaks around $o ~ 2.32 and 9\ = 6q + 7r which are points close to double parton 
scattering kinematics. 

information of the numerator and is therefore capable of evaluating the amplitude 
including all rational terms. The tensor reduction library Golem95 has been extended 
with integrals with massive loop propagators; a new version extending the integrals 
further to the case of complex masses will be released soon. We have presented a 
method for the numerical reconstruction of tensor coefficients from the numerator of the 
integrand of a one-loop Feynman diagram, which can be used to combine the advantages 
of unitarity based and conventional tensor reduction methods. These different tensor 
reduction strategies are combined in the one-loop matrix element generator golem-2 . 0, 
which has been used for the calculation of the QCD corrections of the process pp — > 
bbbb + X. Phenomeno logical results for this process are to be expected soon. 
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